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We calculate predictions from wide class of 'active' models of cosmic structure formation which 
allows us to scan the space of possible defect models. We calculate the linear cold dark matter power 
spectrum and Cosmic Microwave Background (CMB) anisotropies over all observable scales using a 
full linear Einstein-Boltzmann code. Our main result, which has already been reported, points to a 
serious problem reconciling the observed amplitude of the large-scale galaxy distribution with the 
- COBE normalization. Here, we describe our methods and results in detail. The problem is present 

I for a wide range of defect parameters, which can be used to represent potential differences among 

. defect models, as well as possible systematic numerical errors. We explicitly examine the impact of 

varying the defect model parameters and we show how the results substantiate these conclusions. 
The standard scaling defect models are in serious conflict with the current data, and we show 
how attempts to resolve the problem by considering non-scaling defects or modified stress-energy 
components would require radical departures from what has become the standard picture. 
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> ■ I. INTRODUCTION 



Topological defects are an almost generic phenomena in nature and have been already detected in a number of 



laboratory systems (see, for example, ref. |1 



where symmetry breaking phase transitions take place. Probably 



the most exciting possibility, however, is that they are formed during spontaneous symmetry breaking at a phase 



(N 



t — . transition in the early universe l^^l, since they could act as the primordial seeds for galaxy formation; the most 
0^ ' plausible models being the so called cosmic string |^,|| and cosmic texture |7| theories. 

(— I [ Theories for galaxy formation can be described as either 'passive' or 'active' |^. In passive theories, such as those 
Q^' predicted by the inflationary paradigm, all the perturbations are set up effectively as (super-horizon) initial conditions 
at very early times, which then evolve under a deterministic linear evolution, until very late when non-linear processes 
take over on the very smallest scales. By contrast, perturbations are created on all scales at all times in active models 
making predictions much more difficult to calculate. Typically, one has to deal with the fundamental non-linearity of 
, the source over a large dynamic range — approximately 25 orders of magnitude — from the defect formation to the 
present. 

The last few months have seen dramatic progress in pinning down the predictions from defect models of cosmic 
structure formation in what we shall describe as the standard scenario, that is defect motivated stress-energy compo- 



X, . . 

5h nents with an assumption of perfect scaling from formation to the present day. Three groups pHTll have performed 
. . . ' calculations which integrate the linear Einstein-Boltzmann equations using the latest technology for different models 
of the defect stress-energy two-point functions to produce predictions of power spectra for the cold dark matter (CDM) 
density field and the Cosmic Microwave Background (CMB) anisotropies. This article gives a detailed presentation 
of the methods and results of ref. . 

There are two traditional approaches to the study of defect dynamics. Some authors have used large-scale simula- 
tions to provide the sources for their CMB and structure formation calculations ||l^-|lq], while others have developed 
analytic models which attempt to describe the statistical properties of the defects ]lq-Q8[ . Even using the latest 
technology, modelling the source using simulations is severely constrained by dynamic range. In this work, we will use 
a model based approach to calculate the two-point correlation functions, which act as sources for a state-of-the-art 
linear Einstein-Boltzmann solver. This has a number of advantages and also disadvantages when compared with the 
complementary simulation based approach; the real strength being that one is not constrained to a particular defect 
based scenario, allowing one to explore all possible scenarios and understand the robustness of any claims that one 
might make. The downside is that one must take care to construct the source stress-energy, which has many possible 
degrees of freedom, in a way which is at least physically plausible. This is generally done by comparison to some kind 
of simulation. 

The recent work by Pen et al Q and our collaboration |^ , has exposed an apparently deep conflict between COBE 
normalized defect models and the observed galaxy correlations on large scales. The conflict was quantifled in ref. 



in terms of 5ioo, the bias between the dark matter and baryon distributions on scales of 100/i~"'^Mpc required to 
match the COBE normahzed defect models with the galaxy data. Although some evidence for such a problem has 
existed for some time, uncertainties about whether the computations had sufficient dynamic range meant the precise 
quantitative details of the problem were not fully understood. Most previous work on this problem relied on separate 
calculations for the large-angle CMB, which is normalized to COBE, and the linear matter power spectrum, using 
analytic expressions to relate the two, via, for example, the mass per unit length in the case of cosmic strings. Our 
calculations (as well as those of Pen et aZ |^) do not use any such extrapolation, with the perturbations in the matter 
(dark and baryonic) and the photon (plus neutrinos) distributions being calculated in a self consistent way across all 
observable scales. Hence, a single normalization to COBE gives the normalized linear matter power spectrum. 

We should note that the scale of 100/i~^Mpc was chosen for three reasons. Firstly, 100/i~"'^Mpc is sufficiently large 
that non-linear effects are not expected to affect the power spectrum. Secondly, the discrepancy in the power spectra 
at this scale is at its most extreme. And finally, it is unlikely that the power spectrum on these scales can be affected 
by changes in the cosmogony, for example, the introduction of hot dark matter (HDM). 

The price of our solution to the dynamic range problem is that (1) we can only calculate the power spectrum of the 
matter and CMB since we have only included the two-point functions of the source — no non-Gaussian effects 
are included — and (2) the results depend on the validity of the simple scaling picture over many orders of magnitude 
of cosmic expansion. Even though there is substantial support for this assumption, both from numerical simulations 
and analytic modelling pO|- p^ , p^| -fi8[ , there are also reasons to think that it may not be perfect. One of the important 
components of our work is an investigation of possible deviations from scaling. Our approach to modelling the defect 
two-point functions has also allowed us to explore other variations, besides deviations from scaling. These variations 
represent possible differences between defect models and possible systematic errors in numerical defect simulations. 
We have found that solving the 5ioo problem requires extreme departures from the standard picture. (Interestingly, 
we are learning that defect networks in ^Imatter < 1 cosmologies may exhibit interesting levels non-scaling behaviour 

Ha-) 

In the next section, we discuss some aspects of the linearized Einstein-Boltzmann solver, in particular the inclusion of 
source stress-energy, the Einstein-Boltzmann equations for vector perturbations and ways of calculating the ensemble 
average for incoherent perturbations. Section III presents in detail our modelling of the source two-point functions, 
by reference to a specific string motivated model. The results for the standard scaling model are presented along with 
quantification of the 6ioo problem. We also show that simple modifications to the model and cosmological parameters, 
such as the Hubble constant and baryon density, have little impact on the &ioo problem. In Section IV we discuss 
possible deviations from the standard scaling assumption and section V explores further modifications to the model, 
which might lead to an improvement in the comparison with the data. We find that it is possible to get better 
agreement, although only the most extreme modifications come close on 100/i~^Mpc scales. The aim of this paper is 
to present a pedagogical exposition of our work, so that the expert can reproduce and interpret our results. In the 
final section, we discuss its relation to that of others and point to directions of future investigation. 

II. LINEAR EINSTEIN-BOLTZMANN SOLVER 
A. CMBFAST 

In order to calculate CMB anisotropics, one must solve the linearized Einstein-Boltzmann equations. Recent years 
have seen techniques developed to solve these equations to very high precision (~ 1%) for the standard adiabatic 
models based on inflation j2^, culminating in the public release of CMBFAST which can perform this task 
in under a minute on a modern workstation. The standard Boltzmann method involves evolving over 3000 highly 
oscillatory linear ODEs from some time deep in the radiation era to the present day, which can take many hours. 
The line of sight method used in this code reduces the time drastically by splitting the prohibitively oscillatory 
geometric effects from the dynamical processes due to, for example, the Doppler effect and potential evolution. This 
reduces the number of ODEs down to about 30, but adds an integration along the line of sight. 

It is usual to express temperature anisotropics in terms of a decomposition into spherical harmonics, 

^rp OO I 

— [9A)^Y. E ^ImYlrrXOA)- (1) 
i=0 m=~l 

For a Gaussian random field, such as those generated by inflation, the statistics of anisotropics are entirely specified 
by the angular power spectrum C; = where the angled brackets denote an ensemble average. 
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For the moment, let us assume that we only require the anisotropies for a simple inflationary model which creates 
no appreciable vector and tensor fluctuations. In this case, the angular power spectrum is given by 

Ci^^l edk{Af{k,ro)Af*{k,ro)) , (2) 

where the photon distribution function Af (fc, tq) at the present day r tq is given by the line of sight integration, 

Af (fc, to) = / dr {SUk, r)jf [k{ro - r)] + S^k, r)j7° [k{To - r)] + S^k, r)jf ° [k{To - r)]) , (3) 
Jo 

and S'^(fc,T), S^{k,T), S^{k,T) can be deduced from ref. with 

jf'ix)=Mx)^ f°ix)=jlix), jrix) = ^{3jlix)+Ux)) . (4) 

For the coherent limit, implicit in phase focused |^ inflationary models, one can perform the ensemble average by 
simply replacing (Af (fc, tq) Af * (fc, tq)) = |Af(fc,To)|2. 

In the rest of this section, we will describe how this approach can be modified to include active sources such as 
cosmic defects. Firstly, we show how simple coherent scalar sources may be added. Then we discuss the inclusion of 
the vector and tensor sources, almost generic in any active model. Finally, we show how one may perform the ensemble 
average when the source is not coherent and discuss the various implications of decoherence. We have already noted 
that for non-Gaussian sources, such as the topological defects under consideration here, the angular power spectrum 
does not entirely specify the nature of the of anisotropies. However, most realistic models are thought to lead to 
only mildly non-Gaussian anisotropies through the Central Limit Theorem for the superposition of non-Gaussian 
probability distributions. Hence, it should be a useful discriminant between different models for structure formation. 
We will discuss the efficacy of using power spectra to distinguish between different models of structure formation in 
the conclusions. 



B. Coherent active scalar sources 



As a first step, therefore, let us introduce an independent covariantly conserved component of stress-energy 0^^ 
into the Einstein equations, 

G^. = Stt [T^, + e^,] , (5) 

where T^i^ is the stress-energy of CDM, baryons, photons and neutrinos present in a particular cosmogony. This 
extra component, usually assumed to be 'stiff', that is, unaffected by gravity at first order, adds a forcing term to the 
Einstein-Boltzmann equations which represents the active sources. 

For the moment let us assume that the Fourier transform of the stress-energy can be decomposed as 0oi = Qoki 
and 

6,, - y^jQ + (^kkj - , (6) 

where ki is a unit vector in Fourier space, Qo is the velocity field, Q is the isotropic pressure, or three times the 
pressure, and is the anisotropic stress. The conservation equations for this decomposition are 

eoo + - (eoo + e) - = , 

a 

ez3 + -ez5 + ifc2(e + 2e^) = 0. (?) 

a 6 

Hence, in order to incorporate this source stress-energy, one must add to the ODE solver, modify the linearized 
Einstein equations to include the forcing terms and specify two quantities from Qqq, Qd, 'S) and . 

The initial conditions for the Einstein-Boltzmann equations must also be modified. The idea is that one sets up 
initial conditions for a pure growing mode deep in the radiation era. In order to enforce causality, one requires 
that components of the pseudo-stress-energy tensor r^i/ be zero which creates a balance between the initial 
metric perturbations, the defect stress-energy and the matter perturbations. However, there is a residual degree of 
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FIG. 1. The angular power spectrum for the simple coherent sources (a, solid line) and (b,dotted line) normalized to COBE. 
Note no vector or tensor components are calculated here. 



freedom which allows us to just simply set everything to zero. Physically, this implies that the initial conditions are 
unimportant relative to the actual sources themselves, which is implicit in the distinction between passive and active 
sources. 

The most common defect based models are thought to scale for most of the history of the universe and so for 
the moment at least, we specialize our discussion to scaling sources. This requires that T^/^9oo(fc, t) — Fiikr), 
T^/^e{k,T) = F2{kT), T^/^eD{k,T) = F3(fcr) and T^/'^e^{k,T) = Fi{kT), with the functions F,{x) having well 
defined power series expansions around x = 0. Moreover, further constraints can be placed on the leading order 
behaviour of these functions by causality and analyticity In particular, this implies that, 

(ee*) = vKfco + ocfc^), (Oi^ei,) = Ffc'* + o(fc6), (8) 

with all the other correlators and cross-correlators being deduced in a similar way or by assuming stress-energy 
conservation. In the coherent hmit (Ooo©'^*) = (0oo0oo)^^^(©'^©'^*)^^^i which using (||) imphes that F = and in 
fact {Q^Q^*) ^ k"^. This very specialized limit leads to some slight subtleties which will not in general be present for 
active sources. This will be discussed in the section on incoherent sources. 
One simple choice which is consistent with the coherent limit is to define pO[ | 

, ^ 3a sin Akr , 6 / sin Bkr sin Cfcr\ 



-1/2 AkT ' ri/2 B^-C^ \ BkT CkT 

where A, B, C, a and f3 are constants. The angular power spectra for this source are presented in Fig. |]for (a) a = 1, 
(3 = and A = 1 (h) a = 1, f3 = 1, A = 1, B = 1 and C = 0.5. Taking into account the arbitrary normalization used 
in ref. the results seem to be identical. 

At a very technical level, various modifications to internal systems parameters of CMBFAST were also required. 
The source function for Af (fc,T) is dependent on three variables and the code sets up a discrete three-dimensional 
array containing these values with spacings A^, AA: and At, before performing the line of sight integration over I 
values. Obviously, the most accurate results are obtained by using the smallest possible values of Al, Ak and At; the 
cost being an increase in the CPU time. We systematically reduced all three spacings from the values in the original 
version and found that the results converged when Ak was reduced by about a factor of two, with Al and At being 
left the same. We also varied the number of equations solved by the ODE integrator and found that there was no 
discernible improvement from any increase. We concluded that the line of sight integration method is also useful to 
calculate the predictions of active models for structure formation. 



C. Generalization to include vector and tensor modes 



The split of the energy momentum tensor (|6|) is not the most general since it includes only scalar sources. A more 
general split is 

= Qoh + Qo^ , = ^,1^ + (^k " j-^y ) 6^ + [k,Bj + fc, ) + 9^ , (10) 
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which includes vector {Qq^,QY) ^'^d tensor {Qjj) sources such that 

kQl = keY = kQl = kjQl = el^o. (11) 

Without loss of generality, one can fix the direction in Fourier space and choosing k = (0, 0, 1) gives 
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(12) 



with Qq^ = = ©3i = ©?3 = for i = 1,2, 3 and + Q22 — 0- There appear to be four independent vector 
components (e^, e^, eg'i, eQ2) and two tensor components {Qii, ef^^)- However, two of the vector components are 
related to the other two by stress-energy conservation, and hence there are also just two independent components. 

These vector and tensor sources can be identified with their respective contributions to the angular power spectrum 
Ci = Cf + + Cf which are defined in analogy to @ by 

Cl^ - [ k'dk{Al{k,To)Al*{k,To)), (13) 



where the index / runs over I = S,V,T corresponding to the scalar, vector, tensor contributions. The scalar photon 
distribution function is given in (ph, with the vector and tensor versions defined by 



(fc. To) = / dr {V^ik, t)jI' [fc(To - r)] + F|(fc, r)jf ^ [fc(To - r)]) , 

Af{k, to) = r dr {T^{k, r)f' [k{ro - r)]) , (14) 
Jo 

where Vji, Vy, are the source functions for the vector and tensor perturbations and 



The Boltzmann equations for the scalar and tensor perturbations are well studied both analytically and numeri- 
cally, but the same is not true for vector perturbations. Here, we include the Boltzmann equations for the vector 
perturbations, along with a discussion of their salient features. For a more in depth study of scalar, vector and tensor 
perturbations, the reader is referred to ref. ]27t which also includes all the other relevant Boltzmann equations. 

We split up the photon distribution function for the vector sources into its angular multipole moments 

00 

A^(fc,To,M) = EH)'(2; + l)P/(Ai)Ar(A:,To), (16) 
where n = cos 9 is the angular variable. With this decomposition the Boltzmann equations become 

for / > 0, where k is the differential cross-section due to Thomson scattering. 

Si = Ik {Vb + v) , 82 = = ^k (a^ - V6EY) , (18) 
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and S"; = for ^ > 2. In the above expression, is the quadrapole of the electric component of the photon 
polarization distribution and V is the vector metric perturbation which satisfies 

where 0^ is the vector source, either Q\ or 0^, is the average photon density, pn is equivalent quantity for 
neutrinos and N2 is the quadrupole component of the neutrino distribution. Vb is the vector (vortical) component 
of the baryon velocity which satisfies 

= y - ^ (Vb - 1/) + ^ (3A^^ - Vb) , (20) 

where R = 3pb/4:pj and pb is the average baryon density. As already discussed the standard approach is to solve these 
linear differential equations equations for I < Ima.^, plus the equivalent equations for the polarization distribution 
functions, which we have not included here for brevity. The line of sight integration method which we use only 
requires Imax ~ 7, hence reducing the amount of CPU time required. However, we need one final ingredient, the 
source functions for vector perturbations, which are given by 

Vr^ = ke-'^ {Vb - V) , 

V^ = ke-'^P^ + -^e-'^kV , (21) 
V3 



k(t) = / K(T')rfr', (22) 



where 



is the optical depth due to Thomson scattering. We should note that exactly analogous expressions exist for both the 
electric and magnetic components of the polarization. These are also presented in ref. [ pTf . 

The important thing to notice is that vector perturbations are not created at any significant level in the absence of 
a source. This can be seen from the Einstein equation: if Qy = then V oc since the effects of and are 
negligible. Hence, they are not present in inflationary models, but are generic in any active model. 



D. Incoherent sources 



In the previous few sections we have studied how the linearized Einstein-Boltzmann equations can be modified to 
include active sources which create scalar, vector and tensor anisotropics and we have shown how it is possible to relate 
to the measured quantities, such as the angular power spectrum, when the sources are coherent. Such sources are, 
however, very unnatural, generically relying on assumptions such as pure spherical symmetry to maintain coherence. 
Therefore, we are forced to consider methods for calculating the ensemble average for theories which include incoherent 
sources. 

This subject was studied in a series of works which first discussed the equally unnatural, but probably more 
realistic, totally incoherent limit a nd then multi-parameter models which allowed for a gradual relaxation from 

total coherence to total decoherence p'^. One key advantage that the current work has over, for example, ref. js^ is 
that once a given source model is chosen there is no uncertainty in relative normalization of the anisotropy created by 
acoustic oscillations at the time of last scattering and that created much later by the ISW effect. The methods used 
in the earlier work suffered from large uncertainties in this aspect of the calculation. We shall see that decoherence 
is actually prevalent in most active models and the contribution to the anisotropy from the surface of last scattering 
does not, in general, give rise to noticeable acoustic type signatures in the angular power spectrum. We also find a 
suppression of the anisotropy created at last scattering with respect to the ISW effect, not anticipated in refs. |^,^,^. 

The main problem in dealing with decoherence is that the unequal time correlators (UETC) do not factorize, that 

is, 

(P(fc,ri)P*(fc,r2)) ^ {P{k,T,)P*{k,T,)y/^P{k,T2)P*{k,T2)y^\ (23) 

for some arbitrary quantity P{k, r) and ri 7^ T2. There are two methods which one can use to overcome this problem. 
The first |^,|^, used in ref. is to treat the UETC, evaluated at the discrete times used in the linear Einstein- 
Boltzmann solver, as a matrix which is symmetric and hence diagonalizable. This diagonalization yields a change of 
basis and the source can then be written as a sum of coherent sources. 
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(P(fc, Tl)P*(fc, T2)) = TlX (fc, T2) , 

i 



(24) 



where the Ai are the eigenvectors, ordered such that Ai > A2 > .. > A„, and Vi{k, r) are the orthogonal, coherent basis 
functions. Since everything is hnear, one can use each of these basis functions, as a source in the modified version of 
CMBFAST and the sum of the individual angular power spectra yields the total. Although each of these coherent 
sources is has a degree of acausality, it appears to give good convergence from using only the ten largest eigenvalues 
, which indicates that the full calculation is well behaved. 

The second method ||], which is used in ref. and the current work, is to not work directly with the UETC. 
Instead, we create an ensemble of source histories which has the same two-point correlation statistics as the required 
UETC. The exact process for doing this is discussed in the next section, but once it is done, the ensemble average 
can be estimated by averaging the angular power spectra from many individual source histories. Since our ensemble 
is finite, one can also calculate its standard deviation to give an idea of how accurate the calculated average is. The 
results presented in this work used either 100 or 400 source histories to give small statistical errors. However, it was 
possible to gain a qualitative feel from as little as 40 source histories. 

This method is clearly more computationally intensive than the eigenvalue decomposition requiring CMBFAST 
to be run over 100 times as opposed to about 10. We improved the turnaround speed by doing the calculations in 
parallel and 100 realizations took about an hour on IGxRlOOOO processors of the new SGI Origin 2000 owned the UK 
Computational Cosmology Consortium. Setting aside the apparent computational inefficiency|^, we believe that this 
method is more physically transparent since each of the source histories will be causal and at the very least provides 
a useful check on the diagonalization method. 

For these incoherent source histories the nature of the anisotropic stress sources will be seen to be important on 
super-horizon scale, since it is implicitly linked to the vector and tensor sources. As mentioned earlier, in the special 
case of the coherent limit (O'^O'^*) ~ fc*, which is not true in general. Using similar arguments, it was shown in ref. 
!§, that 

(|e^P):(|erP):(|e^f)^3:2:4, (25) 

where all indices are summed. The reason for this is that each of the components is linked via the anisotropic part of 
the space-space component of the source stress-energy tensor. Using a simple model, it was deduced that 

Cf : CY ■■ Cf = 1.46 : 1 : 0.29 , (26) 

at around ^ = 10. While it is true that, the vector and tensor contributions to the angular power spectrum are likely 
to be similar to the scalar contributions on large angular scales, there is no general formula or constraint for this ratio 
and indeed, the model we present in the subsequent sections will be seen to have a larger scalar component than that 
of refs. P,p9t. We should note that the incoherent case, where {Q^Q'^) '--^ fc°, the formation of anisotropics along 
the line of sight, that is, the ISW effect, will be influenced by super-horizon correlations, in a way not possible in a 
coherent model. 



III. THE STANDARD SCALING MODEL 
A. Modeling the source histories 

Here, we present the model for defect two-point functions based on a description of scaling cosmic strings. First, 
we set out the motivation for the model, then we give the mathematical details, explaining each of the parameters. 
Finally, we show some two-point functions calculated within the framework of the model, demonstrating an acceptable 
level of agreement between these two-point functions and those measured in simulations. 

The starting point for our model is ref. p5[ . In this work, measurements of the string two-point correlation 
functions in Minkowski space simulations of network evolution were made and a strikingly simple analytic model was 
put forward, capable of reproducing the important features of these two-point functions with good accuracy. The 
basic assumption of this model is that a string network can be represented as a collection of randomly oriented straight 



^Obviously there will be some overhead in doing the eigenvalue decomposition for any particular UETC and also in creating 
the UETC from simulations. Therefore, we are not quite comparing like with like. 
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segments, each of length ^t, where t is the physical time and ^ is a constant parameterizing the coherence length of 
the string. To model the motion of the strings, each of these segments is assigned a randomly oriented velocity whose 
magnitude is chosen from a Gaussian distribution with zero mean and standard deviation v. 

Under these assumptions, an analytic expression can be derived for the string energy two-point correlation function 
as follows: the Fourier transformed stress energy tensor of a string in Minkowski space is given by 

e^,(k, t) = fijd<j (x^x, - x;x;) e'k.x(.,t) ^ ^27) 

and hence the energy two-point correlator is given by 

(eoo(fc, ti)eSo(fc, t2)) - ( J daida2e^'^-(^(--*i)-^(--*^))) . (28) 

If one now assumes that the quantity (X(cri,ti) — X(cr2,i2)) is Gaussianly distributed, with mean zero and variance 
r, then it follows that 

where a+ — ai + a 2 and cr_ = a\ — ai- To estimate the variance F, one makes use of the idea that on scales smaller 
than the correlation length the string network resembles a collection of straight line segments with velocity ti, which 
implies that 

V[a-MM - {\-v')G''_^v^{t\~tlf , (30) 

for |(T_| < and on scales larger than there are no correlations, so that 

F(CT_,ti,t2) = 00, (31) 

for |cr_| > In this picture, the length and number density of string segments does not change, so that scaling 

behaviour will have to be imposed on the correlators later by hand. Substituting for F in (|2^), performing the 
integration over i7_, and using i J daj\- = VS^^^ (assuming a length of string V^~^ per simulation volume V), one 
finds that 

(e„o(M.)e;„(M.» ^ ^jl^ji^g^f (^|) . (32) 

where erf (x) is the error function, 

erf (a;) = / dx' e-'''\ (33) 



The expression for the energy two-point correlator in (^2|) has the wrong scaling behaviour, since changes in string 
segment length and density have not been incorporated. It is possible introduce the right scaling behaviour into this 
model by hand using 

(eoo(fc,ti)eSo(fc,t2))sca = -=^V^(fc,ii)P(fc,t2)e-^''='(*-*^)' , (34) 



where P{k,t) is times the equal time energy correlator from (32), that is, 

ppik,t) ^ ^t{eoo{k,t)Q*oik,t)) . (35) 

Similar expressions can be derived for the all the other equal and unequal time correlators. 

This model has a number of shortcomings, primarily because certain simplifying assumptions have been made in 
order to make it possible to derive analytic expressions for the two-point functions. By contrast, we need not work 
with simple analytic forms, since we do not work with the unequal time correlators directly. Instead, we use numerical 
techniques to generate histories for the source functions with the correct two-point statistics. This has made it possible 
to improve and extend the model in a number of ways. In particular, we include an improved treatment of causality 
and scaling, and extract a different set of components of the string stress-energy. We briefly sketch these differences, 
before embarking on a detailed mathematical description of our model. 
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• Causality: One problem with the analytic model is that it does not fully respect the constraints imposed by 
causality, which require that there can be no correlations between source components at space-time points whose 
past light cones do not intersect. In particular, by assuming that the lag X((Ti ) — X((T2) between two-points on a 
string segment is Gaussianly distributed, one assigns a non zero value to the probability of correlations existing 
on scales larger than the causal horizon, making the model manifestly acausal. For this reason the oscillations, 
which should appear generically in the two-point functions of causal theories [^3| , are not present in the analytic 
expressions for the correlators. By contrast, our model is causal by design, because we do not assume that this 
lag is Gaussianly distributed. 

• Scaling behaviour: In ref. p5| ], the unequal time correlation functions were simply multiplied by appropriate 
factors by hand in order to enforce the correct scaling behaviour. We extend the model by making the scaling 
form of the unequal time correlators arise in a natural way, as a consequence of the decay of the string segments. 
It will be seen that a different super-horizon form for the unequal time correlators arises as a result of this 
assumption. 

• Choice of stress-energy components: Although we make use of the same basic picture of the string network 
in order to calculate the two-point functions, we extract a different set of components of the string stress- 
energy tensor, calculating the others to maintain stress-energy conservation. We extract the energy Qqq and 
the anisotropic stress , as opposed to Qd, and since we also include tensor and vector contributions in our 
calculation, we also compute vector and tensor source components and 6"^. As discussed in ref. ||], this 
particular choice of scalar components is very natural, as the remaining two components of the stress-energy 
tensor are found to be well behaved on integration of the conservation equations, which is not necessarily true for 
the choice Qqo and Qd- Also, taken in conjunction with 0^ and O"^, these components specify the super-horizon 
perturbations in the most direct manner. 

Having outlined the main differences between our model and that described in ref. ]35| ] , we now proceed to set out 
the specific mathematical details. For a general network of strings in an expanding universe, the stress energy tensor 
has the form 

e^,(x, r)=fij da (eX^X, - e'^X^X:,) ^(3) _ X(a, r)) , (36) 

where fi is the string mass per unit length, X are the coordinates of the string world sheet, parameterized by 
conformal time t and spatial variable cr, dot and prime represent differentiation with respect to r and a respectively, 
£2 = X'V(1 - X2), and S'-^^x.) denotes the D irac delta function in three dimensions. Note that now we are working 
in an expanding universe, each of the string segments will have size ^t. We are interested in histories for the Fourier 
transform of the string stress-energy tensor, which is defined via 

e^^(k,T) = J (i3xe* ''e^^(x,T) = dae'^-^^"'^^ (eX'^X" - e-^X'^X^^j . (37) 

Our conceptual 'string network' consists of a collection of straight line segments, each with an individual label 
m, which 'decay' in a smooth way, completely vanishing by some final time r™. A history for the evolution of the 
complete string stress-energy tensor is then written as a sum of the histories for the stress-energy tensors of the 
individual segments, 

e^.(k,r)=5]e;r,(k,T)r°ff(r,T/,L/)r°"(r,r/,Li,L?) . (38) 

m 

The function is a smooth segment decay function, chosen so that the segment starts to disappear at Lfr'p, 
and has disappeared completely at r™, with the additional features that the stress-energy and its time derivative 
are continuous at ^/t™ and r™, which are necessary in order for the ODE solver to function properly. With these 
properties in mind, we chose the following form for T°^: 

{1 . . .T < LfTf 

l + l(a;3_3x) ...LfTf<T<Tf (39) 
■■■T>Tf 

where 
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x = 2^^f^-l. (40) 

Similarly, r°" is a smooth segment appearance function, with almost identical properties to except that it 
represents the smooth turning on of the segment at early times. By analogy to , we chose 

r . . . r < LjTf 

T°^\r,rj,LlL^)^\ I + liSy - y^) . . . Ljrf < r < I^y (41) 

where 

This function is only included for computational efBciency, since it is possible to ignore any particular string segment 
at times earlier than Ljrf, provided Lj and Lj are sufficiently small. We checked this by varying the values of Lj 
and and found that these variations make very little difference to the total stress-energy provided the values are 
small enough. This is because at any time the stress-energy tensor is dominated by strings whose decay times lie in 
the near future. We choose values for Lj and which are small enough that results are not changed by any further 
decrease. 

During the generation of a particular string history, it is not practical to keep track of every piece of string in our 
conceptual simulation volume. This is because the number density of strings n{T) scales like t~^, so that to have of 
order one string segment remaining by the final simulation time tq , the number of strings we would need to follow from 
the initial simulation time would be of order (ro/Ti)"^. In the case of a mode tracked from before radiation-matter 
equality to the present day, this would require us to follow of order 10^^ strings. Instead, since the 

Nf^V[niTf)~n{Tf + dTf)] (43) 

strings decaying between times t/ and r/ -|- drf in our conceptual 'simulation' volume V are randomly located, we can 

1/2 

replace them by a single string, whose amplitude is multiplied by iV^ ; the power of 1/2 coming from the fact that 
random locations in real space correspond to random phases in Fourier space, so that the amplitude of the Fourier 
transform of a number of such segments sums as a random walk for all k =/= 0. 
The equation for a single source history then becomes 

e^.(k,r) =^1/2^ [n(r7-i)-n(r7)]'^'e- (k,r)T°«(r,r/,L/)T™(r,r/,LhL?) . (44) 

m 

For each source history, we use individual string segments, with values of r/ equally spaced on a logarithmic scale 
between and i^maxTo, where i^maxTo must be later than the final simulation time tq in order that all strings inducing 
significant perturbations at time tq are included^. 

The Fourier transform for each individual string segment is given by 



e™ (k, r) = M dae'^-^- (eXi^Kn - ^"'^ t^'l) , (45) 



where is the length of the string segment at time r and Xm((T, r) are the coordinates of the string world sheet, 
given by 

Xm = + O-X'm + WmTX^ ■ (46) 

For each string segment, Xm is a random location (in practice, we generate k.Xm as a random number between and 
27r), while X'm and Xm are randomly oriented perpendicular unit vectors, such that. 



^The effective total number of strings at any time is given by J2 [^(^7^^) ~ ^i'^T)] ('^^ '^f^^f) Tf,L\,Li) , and the 

normalization of n(r) is chosen to ensure that this quantity is equal to (^t)~^. 
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IX'^I = |X^| = 1, 
0. 



(47) 
(48) 



The string velocity Vm is a random number chosen from a Gaussian distribution with mean zero and standard deviation 
w, truncated to prevent \vm\ > 1- 

Performing the integration over cr, and taking only the real part, we find that 



^rsinc (kX'^(_T/2^ ^cos(k.Xin) cos(kX3Vnir) + sin(k.Xm) sin(kX3V,nT' 



vfjdj - (1 - vDx'a' 



ft'' 



(49) 
(50) 



where sinc(x) = sin(a;)/x and the subscripts refer to the individual spatial components. For conciseness, we have now 
dropped the subscript m on X here, and in the following equations. 

As already noted there are two independent vector and tensor components of the stress-energy, which are sourcing 
the perturbations. However, each of these components will have the same two-point correlation statistics and hence 
we need only evolve one of each and multiply by the appropriate normalization once the power spectra are calculated. 
The anisotropic stress, vector and tensor components are given in terms of the spatial stresses ft^ by 

(51) 
(52) 
(53) 



8^ 


= (2633 


-9ii - 922)/2, 




= er = 


©13, 




= 9^2 = 


9l2. 



For each individual string segment, we find that 



9- = i 



9^ = 
9^ = 



3X3X3-1)) -(l-«,l) (3X^X^-1) 
©00 , 

©00- 



9, 



00 ■ 



vlX.Xs - (1 - vl)X[Xl, 
vlxdi - (1 - vl)X[X'^ 



(54) 
(55) 
(56) 



Integrating over the random orientation vectors, we find that for a single string the super-horizon ratios are in 
agreement with ( |25| ) and since the total stress-energy tensor for the string network is just a sum over the contributions 
from the individual segments, we find that the super-horizon forms of the total stress-energy are also in this ratio. 
However, we have already noted that the ratio of 9^^ to ©00 is not constrained in a similar way, and is likely to be 
highly model dependent. For our model, we find that on super-horizon scales. 



I© 



S|2\ 



l©G 



15w* : 10, 



(57) 



if we assume that the velocities are Gaussianly distributed, rather than the truncated Gaussian which we use in 
practice. This limit, which has been used to make the problem analytically tractable, will be realized for small v. 

Having worked out the energy and anisotropic stress, the remaining scalar components follow by stress-energy 
conservation. By rearranging these equations, we find the following differential equation for ©£> in terms of ©00 and 
©5 



©i 



-2«©, 
a 



fc2 

3 



(^(e„-e„.) 



© 



00 



2©* 



0, 



while © is just 



= 



©D - 9oo I - 9 



(58) 



(59) 



In practice, we use the techniques described earlier to compute histories for the components 9oo, 9'^, ©^ and ©-^. 
Values of each component are stored for a set of times which are closely enough spaced that a linear interpolation 
scheme can accurately reproduce the full history for the function and its derivative. These interpolated functions are 
then used as a set of driving terms to the ODE solver in CMBFAST |2^. In order to increase speed, the evolution of 
©D is only carried out for times satisfying kr < Xmax, where the parameter a;max is chosen to be large enough that 
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further increases do not affect the results, and for later times &d is set to zero. This can be done because in all the 
models we consider here perturbations in Qo are suppressed on scales much smaller than the horizon scale. 

At this point we comment on the way in which stress-energy conservation and compensation are treated in our 
model. In constructing forms for Sqo and we have only been thinking about the behaviour of the long string, and 
not about the behaviour of the loops and gravitational waves into which the long string decays. We have ensured that 
stress-energy conservation is satisfied by only computing two scalar components and using the conservation equations 
to work out the other two. 

One way to treat the loops and gravity waves explicitly is to consider the source G^j/ to be the sum of two 
components, a long string component L^j, and a second fluid component S"^!^. We then model the rate at which 
energy and momentum are being dumped from the long string into the second fluid, which in this case is loops and 
gravity waves, by introducing two functions go and gu, with L^^ satisfying 

d 

Loo H — (-^00 + L) — Lb = —.go , 
a 

Ld + 2-Ld + ^{L + 2L^) = -go , (60) 
a 3 

and Si_ii, satisfying 

Soo H — {Soo + S) — Sd = 5o , 
a 

Sd + 2-Sd + ^{S + 2S^) - go ■ (61) 
a 6 

Given a model for LqO) Ld, L and , such as the one described above, plus an equation of state for the second fluid, 
we can then compute the loop production functions and all the components of S^i, and hence, we have the total, 
8^1,. In the model described here, we have effectively done this by setting 5*00 = = Q. Although this choice does 
not correspond to a particular, identifiable fluid, we have found in studies of the CMB anisotropics created at the 
surface of last scattering, that it gives results which are very similar to physical models for 5'^i/, such as free-streaming 
massless particles. In particular, the main conclusions of this paper and ref. [|l0| will unchanged. However, more 
detailed Modeling will be required if accurate predictions are required. The results of an in depth study of this issue 
will be presented elsewhere ||36|| . 

We now present a sample of two-point functions calculated using these techniques. In the left hand graph of Fig. 
^ we show equal time two-point functions for Bqo and Qdi together with fitting functions for the same two-point 
functions as measured in the simulations of ref. ]3^ ]. The noisy curves are those computed in our model, using 8000 
realizations, while the smooth truncated curves are those of ref. js^. The two-point functions for in this graph 
are obtained by integrating (|5^) for each history. In order to make sensible comparisons between our expanding 
universe calculations and the Minkowski space simulations, we compare our conformal lengths and times with their 
physical lengths and times. Firstly, it should be noted that the dynamic range probed by the simulations is small, 
whereas within the framework of the model the dynamic range can be extended arbitrarily. Within the range probed 
by the simulations, the model appears to give two-point functions in good agreement; the one exception being the 
limiting behaviour of the Qd self correlator. However, it should be noted that the simulations only probe for the 
long string, not the loops and gravitational radiation which the long string spits off. In fact, the fitting function for 
the self correlator has a super-horizon form which is inconsistent with causality and stress-energy conservation, 
since (0ri8J)) ~ rather than fc"*. Our model on the other hand fully respects stress-energy conservation with 
{QdQ*d) ^ fc^, so it is not surprising that there is some level of discrepancy between the limiting forms of the 
functions for this particular component. 

The exact forms of the two-point correlators within our model depends on the choice of string parameters v and 
^. We find that optimal agreement between our two-point functions and those of ref. is obtained when we input 
values of v and ^ which are slightly different to those which are actually measured in the simulations. For Fig. |^ we 
use V = 0.35, ^ = 0.15. In this respect, our model does slightly worse than that reported in ref. |3^], which manages 
to achieve a miraculously good fit to the amplitude and form of the energy equal time cross correlator using exactly 
the values of the parameters v = 0.6 and ^ = 0.15 which were measured in simulations. In spite of this, the limiting 
behaviour of the two-point functions has the correct form in our model and for some choice of the parameters v and 
^ we are able to obtain a good fit to the correlators measured in the simulations. 

In the left hand graph of Fig. || we show equal time two-point functions for 8oo, O'^ and their cross correlator, 
along with 2a error-bars, computed using 8000 realizations. We see that far outside the horizon the cross-correlator is 
relatively noisy, but its behaviour is consistent with a power law of fc^ everywhere except inside the horizon, where it 
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is of order the two self correlators. In fact, it is easy to show analytically within the framework of the model that the 
cross correlator must go like outside the horizon in the limit of a large number of realizations and this behaviour 
clearly manifests itself in the range t = 20 to t — 100. We should note that the noisy behaviour of the cross correlator 
far outside the horizon does not appear to have too large an effect on the matter and CMB power spectra, for which 
the ensemble average has a relatively small variance even for only 40 realizations. 

In Fig. ^ shows the unequal time correlation function for the energy and the corresponding plot from ref. ||35|] . It 
can be seen that the sub-horizon form of the unequal time correlators is similar in both models. However, we see that 
our unequal time correlators have a distinctly different form on super-horizon scales. We quantify this difference by 
using the fimction U(k, ti, T2) defined in terms of the violation of the factorization relation (23) as 



U{k,Ti,T2) 



(P(fc,Ti)F*(fc,T2)) 



(P(fc,Ti)P*(fc,Ti))i/2(P(fc,r2)P*(fc,r2))i/2 



(62) 



for some arbitrary function P{k,T), where ti and T2 are the two times in question, with T2 > ti. In our model, only 
those strings which are present both at ti and T2 can contribute to the cross correlator and hence only those strings 
present at the later time T2 can contribute, implying that (P(fc, Ti)P*(fc, T2)) cx (P(A;, T2)P*(A;, T2)). Hence, we find 



Uik,Ti,T2) CX 



' {P{k,T2)P*{k,T2)) 

(P(fc,ri)P*(fc,ri)) 



(63) 



which outside the horizon gives 



U(k,Ti,T2) (X 



1/2 



(64) 



On the other hand in ref. 
decay, with 



Pq] , the super- horizon fall-off of the unequal time correlators is modelled as an exponential 



C/(fc,Tl,T2) 



— p-(Tl-T2) /t 



(65) 



where the coherence time Tc grows like k~^ outside the horizon. This behaviour gives a good fit on the sub-horizon 
scales which their simulations primarily probe. However, on super-horizon scales the power-law fall-off evident in our 
model must eventually dominate. 

In summary, therefore, we have outlined methods for creating source histories based on a model with two parameters, 
the rms|^ speed of the strings v and the persistence length ^, which are measured in simulations. In doing this we have 
been forced to introduce various 'system' parameters, to allow the problem to be solved in a finite time on a discrete 
system, such as a computer. The value of each of these parameters was chosen, so that further increases or decreases 
toward the continuum value resulted in no change in the two-point functions. In particular, for results presented in 
this paper, we used Xmax = 1000, L} — 0.001, if — 0.01, Ns = 200 and Fmax = 10.0. We have also introduced 
the parameter L^, quantifying the rate at which string segments ar e tur ned off. Unlike the systems parameters, Lf 
clearly has some degree of physical significance. However, in section HI C we demonstrate that the dependence of the 
results on the value of Lf is relatively weak, and we choose to use the value Lf — 0.5 for the rest of our computations. 



B. Power spectra for the standard scaling model 

We define the standard scaling model to be one which uses the above two-point functions with the model parameters 
V = 0.65 and ^ — 0.3, measured in expanding universe simulations Q and an assumption of perfect scaling from defect 
formation to the present day. Also, we must specify a particular cosmogony and we do this by analogy to what 



^As mentioned earlier the distribution of strings has been truncated to prevent strings moving faster than the speed of light. 
This prevents v from being exactly the rms value, the difference from the rms value being minimized for small v. 

* Although, note the earlier comment, that we find better agreement with the two-point functions measured in fiat space 
simulations for slightly different values of u an ^ when we use our causal, stress-energy conserving model. We have decided to 
use the calculated values from expanding universe simulations as our standard since they are likely to be more relevant for our 
model. 
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FIG. 2. The left hand graph shows equal time two-point correlation functions from our standard scaling model: 
T{eoo(fc,r)eSo(fc,T)> (solid line), {QD{k,T)<do{k,T)) /k^ (dashed line), t'^ {Qoo{k,T)<do{k,T)) /k (dotted line). Two-point 
functions measured in Minkowski space string simulations are shown for comparison, using the same line types, but trun- 
cated to illustrate the approximate range probed by the simulations. The right hand graph also shows equal time two-point 
correlation functions from our standard scaling model: T{Ooo(fc, r)02o(fe, t)) (dash-line), T(0^(fc, r)0^*(fc, r)) (dotted-line), 
T{Ooo(fc, 'r)0'^*(fc, r)) (solid-line). On each graph, the x-axis is t and k is 0.1 Mpc~^. 



has become called Standard Cold Dark Matter, that is, a flat background (iltot = 1) spacetime comprising 95% 
coUisionless cold dark matter and 5% baryons {Vic = 0.95, f2& = 0.05), with a Hubble constant at the present day of 
ffo — 50kmsec~^Mpc~^. Fig. ^ shows the resulting power spectra, normalized to COBE, for the CMB and CDM 
(solid lines) compared with the standard adiabatic scenario based on inflation (dot-dashed line) and the published 
data points with error-bars based on the assumption of Gaussianity js^js^ . 

The CMB angular power spectrum appears to have no discernible Doppler peak for two reasons: firstly, there is a 
substantial ISW component to the scalar, vector and tensor anisotropics. The split into the different components is 
illustrated for this model in Fig. || and we see that the scalars are larger than the vectors, with the tensors further 
suppressed relative to the other two. More precisely, we find that the contributions to the angular power spectrum 
are in the ratio, 

Cf -.CY -.Cf = 'i:\: 0.4 , (66) 

at / ~ 10. Although the difference between our models and those presented refs. [p|j29[| are only at the level of a factor 
of two or so, it is still worth noting the discrepancy as a direction for future work. We suggest that this is due to a 
difference in super-horizon ratio of (jO'^P) and (|0ooP), already discussed in an earlier section. 

And secondly, the component of the angular power spectrum created at the surface of last scattering is incoherent, 
with the secondary Doppler peaks being cancelled out by decoherence as suggested in refs. ||,^,^. This leads to a 
further suppression of the amplitude in the ensemble average, relative to the large angular scales, since we are averaging 
high peaks and low troughs. We should note that although the comparison with the published CMB data does not 
appear to be good, the plotted error-bars are only at the level of one sigma and deviations from non- Gaussianity may 
require even larger errorbars, particularly for experiments with small sky coverage. We expect the situation to be 
much clearer when the new CMB data begins to arrive in the very near future. 

However, the situation seems to be much more clear-cut in the case of the CDM power spectrum. Once normalized 
to COBE, the linear power spectrum of the CDM appears to fit the data extremely badly, with the predicted curve 
lying much further outside the observational errorbars than in the case of the CMB angular power spectrum. Again 
the errorbars are based on an assumption of Gaussianity and consideration of a non-Gaussian theory will no doubt 
require us to increase their size, but the level of disagreement is much larger than seems likely in any of the realistic 
scenarios, which are thought to be only mildly non-Gaussian on these scales. If we assume for the moment that we 
can compare the theoretical curves with the data in this very naive way, resolution of the absence of power on any 
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FIG. 3. The left hand graph shows the unequal time correlation function (Ooo(fe, t"i)6oo('^i '''2)) from our standard scaling 
model. The z-axis is (rir2)^''^ (6oo9oo)i the x-axis is Log^Q(kTi) and the y-axis is Logj^Q(kT2). The right hand graph shows the 
same plot created using the model from ref.[30]. 
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FIG. 4. The COBE normalized angular power spectrum of CMB anisotropies, (left-hand graph) and the matter power 
spectrum (right-hand graph) for the standard string model (solid curve). The contributions from the defects during two time 
windows (1300 < z < 100 - Long dashed, 100 < z < 1.6 - short dashed) are also included for comparison. Standard CDM 
(dot-dash curve) and observational data (data points) are included for comparison. 
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FIG. 5. The COBE normalized angular power spectrum of CMB anisotropies, (left-hand graph) and the matter power 
spectrum (right-hand graph) for the standard string model. The total anisotropy (solid line) is illustrated, along with the 
partial contributions from scalar (dotted), vector (short-dash) and tensor(long-dash) components in the case of the CMB 
anisotropies. The errorbars show la statistical uncertainties derived from the finite number of realizations in our calculations. 
Observational data and the prediction for standard CDM (dot-dash curve) are also illustrated in the case of the matter power 
spectrum. The dotted line on the matter spectrum shows the prediction of Ref. as discussed in the text. 
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particular scale requires us to postulate a bias between the CDM and the data. While the idea of a bias between the 
CDM and baryonic matter, probed by the catalogues of galaxies and clusters of galaxies which make up the dataset, 
is not uncommon, large values (> 2) are thought to be unrealistic. 

The bias required on a scale of Rh~^Mpc, where now Hq — lOO/i kmsec~^Mpc~^, can be estimated by calculating 
the fractional matter over-density in a ball of radius i?ft,^^Mpc, 

= 47T f ^k^P{k)\W{kR)\^ , (67) 



where the window function W{x) is given by 

3 

M^fa;) = ^ (sinx — a:cosa;) , (68) 

and comparing it to <T}i calculated from a hypothetical curve which fits the data, that is, bji = (JbJc'r^^ ■ The scale 
R = 8 corresponds to an ~ 1, that is, scales turning non-linear at the present day, and is the most common scale 
on which comparisons are made. For standard CDM as — 1-2 for h = 0.5, while the value favoured by observations 
is as — 0.5, illustrating the celebrated excess of power on small scales for this model. When we perform the same 
calculation for the string model, we get as = 0.31 and hence the bias on these scales is bs ~ 1.5. Such a value is 
around the limit of what is thought to be possible, but given the uncertainties involved is not totally unreasonable. 

This comparison does, however, ignore the obvious woeful absence of power on much larger scales. In order 
to quantify this we choose R = 100 since (1) it is unlikely that scales of 100/i~^Mpc are affected by non-linear 
gravitational evolution, (2) the distribution of galaxies is likely to be more Gaussian than on smaller scales and (3) 
such scales are above the neutrino free-streaming scale and so the introduction of a hot dark matter (HDM) component 
cannot be used to modify the shape of the spectrum. We estimate that the standard string model requires a bias of 
bioo — 5.4 to reconcile it with the data. Since the chances of either the actual Universe |39|j40[] or the physical model 
pl| , ^ having such a bias seem remote, we conclude that the standard string model is in serious conflict with the 
current observations on scales around 100/i~^Mpc. 

The COBE normalization of our standard scaling model also allows us to calculate a value for the dimensionless 
quantity Gfj,, where fj, is the string mass per unit length, and G is Newton's gravitational constant. For our standard 
scaling model, we find that G/i = 2.0 x 10~^, which is very close the value G/x = 1.7 x 10~^ obtained from calculations 
of large angle CMB anisotropics using high resolution local string simulations ||l^ . Although these values are in good 
agreement, we stress that the purpose of our work has not been to compute G/i to high precision and we expect 
that variations in details of the string evolution which we have not attempted account for in our model, such as 
the amount of small scale structure, could give rise to variations in G/i. Instead, we are primarily interested in the 
relative normalization of anisotropics on different scales, in particular between COBE and crioo, which can be obtained 
without knowledge of the absolute value of G/i. In future sections we present a number of variations on our standard 
scaling model, each giving rise to a different value for the string mass per unit length. To emphasize the fact that we 
do not intend to use our model to obtain precise predictions for the absolute value of Gfi in realistic cosmic string 
scenarios, we quote all subsequent values in terms of the ratio /i//is, where /Js is the value obtained for our standard 
scaling model. We should note that large values of n/ns are likely to be excluded by the absence of residuals in timing 
measurements for milli-second pulsars ]4^ ] . 

As one final point, in Fig.^ we include the prediction for the matter power spectrum for strings in a background of 
CDM computed in ref. [Q. The dotted line on the matter graph shows the prediction from the I model described in 
this reference, normalized to give the same value for as as our standard scaling model. We note that the predictions 
are in extremely close agreement, although in the current work the amplitude is no longer arbitrary, but fixed by 
normalization of the associated CMB fluctuations to COBE. 



C. Modifications to model parameters 



In order to test the robustness of the conclusions about the standard model, we have repeated the calculations 
for different values of ^ and v. These variations could either represent fundamentally different types of cosmic 
strings to those modelled in simulations, for example, those with superconducting currents, or any possible systematic 
uncertainties in the measurement of these quantities. Fig. ^ shows the power spectra for different values of ^ with 
V = 0.65, while Fig. |^ shows the power spectra for different values of v with ^ = 0.3, with the values of 6ioo, bs, and 
/<//is summarized in table |. All other parameters are as in the standard scaling model. 

The first thing to note is that, although our choice of values spans the possible parameter space, none of the models 
does significantly better than the standard scaling case as far as the value of bioo is concerned and there is also no 
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Short-dash 
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IC 
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1 no 


= 0.03 
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Qt = 0.1 






Short-dash 


1.87 
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1.00 
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Long-dash 


2.17 


5.94 


1.00 


"Best of all worlds" 


U 


Long- short dash 


0.3 


1.56 


0.03 



TABLE I. Table of biases and values of /i//is for the standard scaling model, and the simplest variations described in this 
paper. Each model is labelled by the figure and line type where it appears. 

discernible Doppler peak in the CMB spectrum. However, at a more microscopic level there are differences which 
follow the trends of previous calculations, giving us confidence that the model is reproducing intuitive results. 

As ^ decreases one is reducing the scale at which the two-point functions turnover from white-noise, with ^ = 2.0 
being close the causal limit Q]. It was predicted in ref. that this would lead to the contribution to the CMB 
anisotropy from the surface of last scattering being peaked at smaller scales, which is indeed what we see in Fig. 
although this feature is slightly masked out by the large ISW component. This can be seen more clearly in the 
CDM power spectrum which turns over at a scale k ~ ^fcoq- One would of course expect a reduction in ^ to result 
in a substantially higher peak, but this is partially counter balanced by an almost equivalent increase in the ISW 
component, the only direct evidence for such an effect being that smaller values of ^ require a smaller value of G/i. 
This change in G/x can be understood from the formula ps ~ m/^^! if the value of ^ decreases the value of [i must also 
decrease to keep the string density and hence the amplitude of the CMB anisotropy the same. 

The lack of dependence on v is less intuitive. One might think that in the limit of w — *■ 0, the network would 
become more coherent, since the strings are not moving, creating perturbations in the same place. However, such a 
limit does not correspond in any way to the standard picture of string evolution. If the strings are moving slowly 
or are stationary, then reconnection will take place only very infrequently and the scaling regime will be difhcult to 
attain. But scaling is implicit in our model, being put in by hand, so even if they are not moving the strings will have 
to decay in some way, which is essentially random, introducing decoherence. The only discernible effect of changing 
w is a weak dependence of the amplitude of the matter power spectrum, believed to be due to the dependence of the 
relative amplitudes of Oqo and G'^ on v already discussed. Also interesting, however, is the apparent independence of 
the turnover of the matter power spectrum, suggesting that the corresponding turnover in the two-point functions is 
also independent of ti. Once again the values of 63, ^100 and y^j y^s are summarized in table |. 

We have also tested the dependence of our results on the parameter L f , representing the rate at which the strings 
are turned off. CMB and matter power spectra for various values of L/ are illustrated in Fig. ||. The results have a 
weak dependence on Lf, with smaller values oi Lf giving slightly better values for 6100 than larger values. We note 
that further decreases in Ly below 0.01 do not change the results further. Since the dependence on L/ is only weak. 
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Description Figure Line type 
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TABLE II. Table of biases and values of p//is for each of the models with deviations from scaling. Each model is labelled 
by the figure and line type where it appears. 
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TABLE TIL Table of biases and values of us for each of the models where our standard source has been further modified. 
Each model is labelled by the figure and line type where it appears. Note that we have not calculated jJ./ jJ.s for coherent models, 
which are not based on string-like two-point functions. 
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1 k/h Mpc 1 

FIG. 6. Varying ^, the string coherence length: We plot the COBE normalized angular power spectrum of CMB anisotropies, 
(left-hand graph) and the matter power spectrum (right-hand graph) for various values of the parameter ^, with v — 0.65. (^ = 2.0 
- long-short-dash line, ^ — 1.0 - dotted line, ^ = 0.3 - short-dash line, ^ — 0.03 ~ long-dash line). Observational data and the 
prediction for standard CDM (dot-dash curve) are included for comparison. 
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FIG. 7. Varying v, the string velocity: We plot the COBE normalized angular power spectrum of CMB anisotropies, 
(left-hand graph) and the matter power spectrum (right-hand graph) for various values of the parameter v, with ^ = 0.3.(i; = 0.0 
- long-short-dash line, v — 0.3 - dotted line, v = 0.65 - short-dash line, v = 0.99 - long-dash line). Observational data and the 
prediction for standard CDM (dot-dash curve) are included for comparison. 
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FIG. 8. Varying Lf, the string 'turn off' parameter: We plot the COBE normalized angular power spectrum of CMB 
anisotropies, (left-hand graph) and the matter power spectrum (right-hand graph) for various values of the parameter Lf, with 
5 = 0.3, V — 0.65. (L/ = 0.01 - long-short-dash line, Lf = 0.1 - dotted line, Lf — 0.3 - solid line, Lf — 0.5 - short-dash 
line, Lf — 0.7 - long-dash line, Lf — 0.9 - dot-short-dash line). Observational data and the prediction for standard CDM 
(dot-long-dash curve) are included for comparison. 



we choose an intermediate value oi Lf — 0.5 for the remaining calculations in this paper (with the exception of fig |Tl| 
in which we illustrate the possible improvement to 6ioo which could result from exploiting all conceivable uncertainties 
in our model). The fact that the results depend minimally on L/ also represents evidence to suggest that the results 
will not depend strongly on the exact way in which the decay of long string is treated. 



D. Modifications to cosmological parameters 



We must also consider the possibility of different cosmogonies, since most cosmological parameters are not con- 
strained to better than a factor two. It is simple to change the Hubble constant and also the relative content of 
baryons and CDM. The resulting spectra are presented in Fig. |^ for h ~ 0.3 to = 0.7 and in Fig. |l^ for ilf, = 0.01 to 
rif, = 0.2, keeping fitot = 1 and using the standard scaling model for the two-point functions. Once again, no model 
significantly improves the value of 6ioo (see table ^ . 

The CMB angular power spectrum is largely unaffected by the changes in cosmological parameters that we have 
tried. However, we do see that the shape of the CDM power spectrum is modified by changes in /i, via the time of 
equal-matter radiation and the well-known shape parameter F k, Wn. This fixes the position of the turnover in the 
power spectrum, with larger values of h leading to a turnover at smaller scales. We also notice slight oscillations in 
the power spectrum for larger values of SI;,. In these models the oscillations that are present in the photon-baryon 
fluid are transferred to the CDM, but in contrast to an inflationary model, they are damped out by the effects of 
decoherence. 

We should comment on the apparent absence of any marked dependence on these cosmological parameters, since the 
CMB spectrum for I > 100 is very strongly dependent on them in the case of the standard adiabatic scenario, and this 
dependence has been suggested as a way of making extremely accurate estimations of many cosmological parameters 
using satellite experiments. To understand this difference, we should remember that the anisotropies created at the 
surface of last scattering by acoustic waves are incoherent, leading an absence of secondary peaks, effectively washing 
out the strong dependence on Qb and Hq. And more importantly, any contribution from the last scattering surface 
seems to be swamped by the ISW effect, which is less sensitive to changes in the cosmogony. 

One modification to the standard scenario which is often used to allow the standard cold dark matter model to 
obtain a better fit to the galaxy data is to introduce a small amount of hot dark matter in the form of neutrinos. 
A similar, procedure would also allow the standard string model to fit the shape of the observed power spectrum 
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FIG. 9. Varying h: We plot the COBE normalized angular power spectrum of CMB anisotropies, (left-hand graph) and the 
matter power spectrum (right-hand graph) for various values of the cosmological expansion rate h, with ^It — 0.05. {h = 0.3 
- long-short-dash line, h = 0.4 - dotted line, h = 0.5 - solid line, h = 0.6 - short-dash line, h = 0.7 - long-dash line). 
Observational data and the prediction for standard CDM (dot-dash curve) are included for comparison. 
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FIG. 10. Varying fit: We plot the COBE normalized angular power spectrum of CMB anisotropies, (left-hand graph) and the 
matter power spectrum (right-hand graph) for various values of the cosmological expansion rate fit,, with h — 0.5. (Sli, = 0.01 
- long-short-dash line, Qt = 0.03 - dotted line, Qt ~ 0.05 - solid line, fib = 0.1 - short-dash line, fib = 0.2 - long-dash line). 
Observational data and the prediction for standard CDM (dot-dash curve) are included for comparison. 
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FIG. 11. The "best of all worlds" models: The COBE normalized angular power spectrum of CMB anisotropies, (left-hand 
graph) and the matter power spectrum (right-hand graph) for the standard scaUng model (dotted line) and two "best of all 
worlds" models. In the first of these models, (long-short dash line), the values of each of the parameters ^, v, h, Qb and Lf 
have been pushed as far as possible in the direction which favours high ctioo, while a standard scaling law is used. The second 
model, (solid line) is identical, except that it also incorporates a fairly plausible deviation from scaling over the matter-radiation 
transition. Observational data and the prediction for standard CDM (dot-dash curve) are included for comparison. 



on scales below the neutrino free streaming scale A 20/i^^Mpc, but we anticipate that this will not be as efficient 
on larger scales, and in particular we expect the introduction of HDM to have little bearing on the 6ioo problem. 
Nonetheless, we plan to investigate the implications of making such a changes in future work. 



E. Summary 

In summary, we have scanned the range of each of the parameters in our model, while maintaining perfect scaling, 
and we have found that none of these simple variations are capable of significantly reducing the deficit of power on 
scales around 100/i~^Mpc. However, some of the parameters do marginally improve the situation, in some cases 
reducing the required &ioo from 5 to ~ 3. One might postulate, therefore, that modifications to all these parameters 
simultaneously might lead to a more substantive amount of power on 100/i~^Mpc scales, and indeed this is the case. 
To illustrate this, we performed a run with ^ = 0.0001, v = 0, h = 0.7, fib = 0.01, and Lf = 0.01, which we describe 
as the 'best of all worlds' model, and the result is presented in Fig. |l^. We see that the situation is improved, but still 
the bias required, &ioo — 1.6, is not unity and we now find a large excess of power on scales around 8h~^Mpc. Also, 
the CMB angular power spectrum appears to be worse fit to the data. While this model serves as a useful caveat to 
our arguments, we believe that pushing the model this far is not realistic within the current understanding of defect 
models. Nonetheless, it may serve as impetus for future model building. 

Except for the caveat described above, the minimal dependence of 6ioo on the wide variations in these parameters 
is already strong evidence to suggest that the 6ioo problem will be a feature of most scaling defect models. In the 
next two sections we further test this idea by examining the results of further modifications and generalizations of 
the standard model, with all deviations being described as perturbations from the standard scaling model. 

We should note that there are two simple variations of the cosmogony which we have ignored in the section on 
cosmological parameters, namely an open universe (fi < 1) or the introduction of cosmological constant (flm+^A = 1). 
We anticipate that these variations will lead to modifications to scaling similar to those described in the next section 
p5|-p8[, and might lead to more acceptable values of 6ioo. An in depth investigation of this problem is the subject of 
ongoing research |23j (see also ref. [Q). 
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IV. VARIATIONS FROM STANDARD SCALING 



A. Motivation and implementation 

In the previous section, we introduced the 5ioo problem for scahng defect models and showed that it is apparently 
robust for a range of different parameters. However, we have also noted that this scaling assumption has only been 
tested using defect simulations with a very small dynamic range. Hence, in the spirit of testing of the standard model, 
we should allow for the possibility of modifications to scaling. In fact, this is the most obvious resolution to the 6100 
problem since the scaling assumption is what relates the contributions from defects on different scales. Deviating 
from scaling would allow us to effectively tilt the power spectrum. A similar approach has become a popular solution 
to the excess power on small scales in standard CDM models based inflation, but in our case we want to create more 
power on large scales. 

This point can be illustrated most effectively by reconsidering Fig. ^ which shows the results for the standard model 
and also the partial contributions from the strings present in during z — 1300 to z = 100 (short dashed curve) and 
between z = 100 and z = 1.6 (long dashed curve). These two curves give us intuition about when the perturbations 
relevant for COBE normalization and 6100 are laid down. For example, if one could create an imbalance between the 
strings present during these two time windows, one could hope to change their relative amplitude and reconcile the 
current data points with COBE normalization. This is possible if one modifies the scaling picture, with most graphic 
illustration being the total removal of the string network around z = 100. 

The first type of deviation we consider is motivated by the mild shift in the behaviour of a string network which 
is observed in simulations under going a radiation-matter transition. Typically, quantities such as the string velocity, 
persistence length, string density, and level of small scale structure are seen to undergo a small change at around 
radiation-matter equality. In general, this shift tends to make strings move slower and be less dense. 

One simple way to implement this step-type transition is to allow the effective energy per unit length to change, 
with the mass per unit length being a factor x larger before the transition than after it. We generate histories for a 
source with normal scaling behaviour and then for each history, we multiply the value of each source component at 
each time r by some factor fR^r). Since we require our source histories to behave smoothly, we implement a smooth 
shift in the value of fn using 

/fl(T) = l + (x-l)T°ff(T,TT,LT), (69) 

where is the same smoothly varying function which was used for turning the string segments off, defined in (|3S|), 
but now LtTt is the starting time of the transition, and tt is the end time. 

We have also tried other ways of implementing a transition, such as introducing a shift in the time dependence of 
the number of strings per unit volume. In the standard case, we have n(r) = {£,tYi^, which we modify by setting 
"iT-iT) = fR{T){S,T)~^ 1 where fa is the smoothly varying transition function given in (|69|). We find that the net results 
from these two ways of implementing the transition are very similar and so in the results section, we concentrate on 
the former, simpler case. 

The second type of deviation we consider is a deviation in the scaling exponent. We implement such a transition 
by altering the dependence of the number of strings per unit volume on time. In the standard scaling picture, there 
is roughly one piece of string per correlation length cubed. Since the correlation length is proportional to the horizon 
size T, we find that the number density of strings as a function of time is, n{T) cx t~'^, which we modify by setting 
n(T) (X T~^^~^'^°'\ with a = being the standard value. Using (|4j) we see that the power spectra of the O^i/'s depend 
on n^r) and hence the time dependence of 800 (which behaves like the square root of the power spectrum) outside 
the horizon is now Oqo = t^(2+") 

B. Illustrative examples 

Figs. |l^, |l^ and |l^ illustrate the results for radiation-matter transition runs (implemented by varying /i) with 
various choices of the parameters x, tt and Lt- The first (Fig. |l^) shows mild transitions, with an amplitude x = 2, 
each lasting for a factor of 10 in conformal time r, that is, Lt = 0.1. Each curve shows results for a different choice 
of final time tt- Initially, we see that as the time of the transition is moved later, the peak in the matter power 
spectrum gets higher and is shifted to larger scales. However, as discussed in the standard scaling section illustrating 
the two time windows, we find that as the time of the transition is increased beyond tt = 1000, the height of the 
peak in the matter power spectrum actually falls again. This is because a very late transition tends to boost the 
perturbations on COBE scales as well as those on scales relevant to the large scale matter power spectrum. It is clear 
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FIG. 12. Varying tt, the end time for the modelled non-scaling radiation-matter transition: We plot the COBE normalized 
angular power spectrum of CMB anisotropics, (left-hand graph) and the matter power spectrum (right-hand graph) for a matter 
radiation transition with various values of the final time tt, with Lt = 0.1 and x = 2.(rT = 100 - long-short-dash line, tt = 400 
- dotted line, tt = 1000 - short-dash line, tt = 5000 - long-dash hne). Observational data and the prediction for standard 
CDM (dot-dash curve) are included for comparison. 




1 10 100 1000 0.01 0.1 



1 k/h Mpc-i 

FIG. 13. Varying x, the factor by which ji^s is assumed to be larger before the radiation-matter transition: We plot the COBE 
normalized angular power spectrum of CMB anisotropics, (left-hand graph) and the matter power spectrum (right-hand graph) 
for a matter radiation transition with various values of the factor x, with tt = 1000 and Lt ~ 0.1. (x = 2 - long-short-dash 
line, X = 5 ~ dotted Une, x = 10 - short-dash line, x = 20 - long-dash line). Observational data and the prediction for standard 
CDM (dot-dash curve) are included for comparison. 
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FIG. 14. Varying Lt, the start time (in units of the end time) of the non-scaling radiation- matter transition modeUed in 
the defect sources: We plot the COBE normalized angular power spectrum of CMB anisotropics, (left-hand graph) and the 
matter power spectrum (right-hand graph) for matter radiation transition models with different values of the length Lt, with 
amplitude x — 5 (tt = 100, Lt = 0.1 - long-short-dash line, tt = 100, Lt ~ 0.8 - dotted line, tt = 1000, Lt = 0.1 - 
short-dash line, tt ~ 1000, Lt ~ 0.8 - long-dash line). Observational data and the prediction for standard CDM (dot-dash 
curve) are included for comparison. 
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FIG. 15. Different types of non-scaling radiation-matter transitions: We plot the COBE normalized angular power spectrum 
of CMB anisotropies, (left-hand graph) and the matter power spectrum (right-hand graph) for two different implementations 
of the matter matter radiation transition, each with x — 5 and L = 0.1. (tt ~ 100, varying n - long-short-dash line, tt = 400, 
varying n - dotted line, tt ~ 100, varying ^ - short-dash line, tt ~ 400, varying fi - long-dash line). Observational data and 
the prediction for standard CDM (dot-dash curve) are included for comparison. 
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FIG. 16. Power law deviations from scaling; We plot the COBE normalized angular power spectrum of CMB anisotropies, 
(left-hand graph) and the matter power spectrum (right-hand graph) for a power law deviation from scaling with various values 
of the parameter a. {a = 0.25 - long-short-dash hue, a — 0.5 - dotted line, a = 0.75 - short-dash line, a = 1.0 - long-dash 
line). Observational data and the prediction for standard CDM (dot-dash curve) are included for comparison. 



that a transition occurring as late as today is equivalent to having no transition at all. None of these reasonable, mild 
transitions (x = 2) significantly improves the &ioo problem. 

In the second figure (Fig. ^3|), we stick to a transition time of tt = 1000, which as we have discussed gives us the 
best chance of introducing a shift between the COBE and large scale matter normalization. Using a transition length 
of Lt = 0.1, we vary the amplitude x of the transition. A value of x = 10 improves things significantly on scales 
of 100/i~^Mpc and x = 20 does better still, but we find that further increases only affect the features of the matter 
spectrum on smaller scales. It is interesting that there is a limiting value for the relative COBE/crioo normalization, 
and that this limiting value happens to fit the observations very well. We note however that the cases which fit the 
large scale matter data require implausibly large values of the transition amplitude Xj for which no precedent has 
been seen in studies of defect evolution. Furthermore, extremely drastic additional alterations would have to be made 
to the model in order to fit the small scale matter and CMB data. 

In the third figure (Fig. 14), we illustrate the dependence on the length Lt of the transition, for two different 



choices of the final time tt, with x fixed to be 5. We see that the length of the transition does not strongly infiuence 
the (Tioo normalization. In fact, we have been unable to find any region of our transition parameter space capable of 
fitting (Tioo, where the choices of parameter values is more plausible than those illustrated in Fig. |l^ 

Fig. |l^ illustrates the comparison between a radiation-matter transition in which /i is varied, with one in which n 
is varied. The first pair of curves shows a transition with Lt = 0.1 and tt = 100. In the case where fi is varied, 
we choose x = 5 (short-dash line) whereas in the case that n is varied, we choose x = 25 (long-short dash line). 
The reason for this difference in amplitudes is that increasing /i by a factor of x increases the power spectrum of the 
perturbations by a factor of x^ , while the same increase in n affects the power spectrum of perturbations by a factor 
of X- We see that the resulting spectra for this pair of models are very similar. The second pair of curves illustrate 
a transition with Lt = 0.1 and tt = 400, for the same choices of x as above. Again, the curves are very similar for 
each of the transition models. Hence, the resulting spectra of perturbations does not seem to be strongly dependent 
on the way in which the transition is implemented. 

Fig. |l^ shows deviations in the scaling exponent of various degrees. We see that models where the deviation from 
scaling is significant enough to bring about a substantial increase in the amount of power in the matter spectrum 
on scales around 100/i^^Mpc are so extreme that they completely miss the small scale matter and CMB data. We 
find that further increases in the scaling exponent beyond a = 1.0 do not affect the uioo normalization, only giving 
rise to significant differences in the resulting power spectra on smaller scales. As in the case of the radiation-matter 
transition, it is interesting to note that there is a limiting value of cioo normalization with increasing alpha, and that 
this limiting value happens to pass through the large scale matter data. 
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As one final point, Fig. ^ shows the results of a relatively mild deviation from scaling (varying /i, x = 2, tt = lOrgq, 
Lt = 0.8) where all of the other standard parameters have been pushed as far as possible in a direction which favours 
a large value for 6100 as in the 'best of all worlds model' = 0.0001, u = 0, ft. = 0.7, — 0.01, and Lf — 0.01) 
discussed earlier. In this case, 6100 is of order one. However, in obtaining a reasonable value for 6100, the resulting 
model totally fails to fit the shape of the matter power spectrum, with an extreme excess of power on smaller scales. 
Although we do not believe this limit will allow the resurrection of the standard defect scaling defect scenario, it does 
present a possible road of attack for the construction of more exotic models which could fit all of the data. 

To summarize, we have presented results for a number of models showing deviations from scaling. None of these 
models is able to fit the large scale matter data without extreme corrections to the standard scaling picture, or forcing 
all other model parameters in the direction which favours large 6100 ■ Even in the cases where the value of 5ioo is 
reasonable, a very considerable amount of further work would need to be carried out in order to make the model fit 
small scale matter and CMB data. 

V. FURTHER MODIFICATIONS TO THE MODEL 

In the previous sections we have discussed possible variations from our standard scaling string model, and we showed 
that only extreme deviations from the standard scaling model can significantly rectify the problem. In section III 
we showed how the 6100 problem is relatively robust to changes in the model parameters v, ^ and Lf (as well as the 
cosmological parameters fib and h) . Although all of the variations we have considered take place within the context of 
our string model, the robustness to these changes already provides evidence to suggest that the results will be similar 
for other types of defect. For instance, the independence of the results on the parameter ^ suggests that results will 
be similar for other types of defect for which the two-point functions cut off at different sub- horizon values. In this 
section we further test this idea by introducing further modifications to the forms of our two-point functions by hand, 
in order to see how extreme these modifications must become before the 6100 situation is significantly improved. 

One change which might alleviate the problem is to alter the relative strength of various components of the stress- 
energy. Recent work |49], which makes use of a coherent approximation to model source behaviour, has suggested 
that defect models with a highly suppressed anisotropic stress might give rise to more acceptable values of the matter 
bias factor b and indeed we have already mentioned that the ratio of O'^ and Ooo is model dependent. In order to 
investigate this possibility in the context of our work, we make a simple modification to our standard scaling model. 
We multiply the energy Oqo by z/, where in models with values of < 1 the significance of the anisotropic stress is 
boosted relative to that of the energy, while in models with > 1, the energy is boosted. We should note that the 
ratio of the anisotropic stress to vector and tensor components is unaltered, as required by isotropy and causality. 

Before presenting results for our full string model, we discuss the effect of these changes in the context of some 
simple coherent models^. In Fig. |l^ we present the CMB and matter power spectra for four such models. The first 
model (which we call CI) has 

"°"-\0 fcT>5, 

es = 0, (70) 



while for the second model (C2) 



600 = , 



For this pair of models, the version with suppressed anisotropic stress has a more acceptable value of 6100 than the 
version with suppressed energy. However, this pair of models is not exactly causal, because the real space two-point 
correlation functions do not exactly vanish outside the horizon. We also present results for the following pair of 
models, which do exactly satisfy this condition. They are (C3), with 



^Since in each of these models, either 0oo or O is zero, su per h orizon constraints on the cross-correlator (Ooo6 ) place no 
constraints on the super horizon behaviour of (see section II B). 
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FIG. 17. Various coherent limits: We plot the COBE normalized angular power spectrum of CMB anisotropies, (left-hand 
graph) and the matter power spectrum (right-hand graph) for a range of coherent models discussed in the text. (Non-causal 
models: CI - long-short-dash line, C2 - dotted line. Causal Models: C3 - short-dash line, C4 - long-dash line). Observational 
data and the prediction for standard CDM (dot-dash curve) are included for comparison. 



sin(fcT) 
kr 



and (C4) with 



es = o, 
eoo = , 



i/2sin(M 

KT 



(72) 



(73) 



The sources in the second pair of models are designed to exhibit sub-horizon decay at a similar value of fcr as 
the first pair of models we have discussed. However, we see that in this explicitly causal case, the model with 
suppressed energy actually gives a more acceptable 6ioo than the model with suppressed anisotropic stress. These 
results represent a surprising contrast to those of the case which is not exactly causal, implying that exact imposition 
of causal constraints is necessary for physically meaningful results. Activity outside the causal horizon presents 
opportunities to seed standard 'adiabatic' perturbations (which will naturally make the 6ioo look better). We suspect 
our slightly acausal model has done this to a remarkable degree. 

In Fig. |l^ we present the results for different values of the energy factor v, modifying our standard string two-point 
functions. As in the causal coherent case above, we see that ctioo is actually better for models with suppressed energy 
rather than suppressed anisotropic stress and in fact no value of v would give rise to a significant improvement in the 
6ioo problem (see table [II). 

Another way in which we have modified our standard model is by imposing a sharper sub-horizon cutoff in the 
source stress-energy. By doing this, we hope to have covered a range of possible defect behaviour including that of 
cosmic textures, whose stress-energy tensor will exhibit a faster sub-horizon fall-off than that of strings, reflecting the 
fact that they have less features on small scales. We implement this particular modification without violating the 
requirements of causality as follows. We introduce a parameter e, and specify that the number of strings with decay 
times Tf satisfying fcry > e is zero, while the number of strin gs w ith decay times satisfying krj < e is unchanged. 
Results for various choices of the parameter e are shown in Fig. 19, We see that 6ioo does not depend strongly on the 
the value of e. 

Finally, we allow modifications to the way in which scaling is enforced in our model. We note that our standard 
scaling model contains two time dependent functions, one being the number of strings per unit volume ^(t), the 
other ^(r) being the mass per unit length of the strings. In the case of standard scahng, n(r) ^ t^^ while /Lj(r) is 
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FIG. 18. Varying v, the enhancement factor for Goo which artificiaUy modifies the ratio between energy and anisotropic 
stress: The COBE normalized angular power spectrum of CMB anisotropies, (left-hand graph) and the matter power spectrum 
(right-hand graph) for standard scaling subject to various values of the parameter v. {v = 1 - solid line - u = Q- long-short-dash 
line, V = —1 - dotted line, p — 0.5 - short-dash line, i/ = oo ~ long-dash line). Observational data and the prediction for 
standard CDM (dot-dash curve) are included for comparison. 
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FIG. 19. Varying the cutoff parameter e: The COBE normalized angular power spectrum of CMB anisotropies, (left-hand 
graph) and the matter power spectrum (right-hand graph) for our standard scaling source, with various values of the parameter 
e, representing the value of kr at which the source is subject to a sharp subhorizon cutoff, (e = 2 - long-short-dash line, e = 5 
- dotted line, e = 10 - short-dash line, e = 20 - long-dash line). Observational data and the prediction for standard CDM 
(dot-dash curve) are included for comparison. 
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FIG. 20. Scaling can be enforced by varying either the string density or the mass per unit length. Here we plot the COBE 
normalized angular power spectrum of CMB anisotropies (left-hand graph), and the matter power spectrum (right-hand graph) 
for a scaling source, with various values of the parameter q, which parameterizes the way scaling is enforced {q = Q corresponds 
to constant mass per unit length). We use v = 0.65 and ^ = 0.3. (g = - long-short-dash line, q = —1 - dotted line, q = —1.25 
- short-dash line, q — —1.5 - long-dash line). Observational data and the prediction for standard CDM (dot-dash curve) are 
included for comparison. 



a constant. We note that an extreme case which would also exhibit scaling behaviour would be a model in which 
the number of strings n was constant, but where /i evolved like t~^/^. Such a case would bear some resemblance to 
the inflationary mimic model of Turok , where there is a constant number of expanding shells, but some feedback 
mechanism is required to modify the shell surface density in order to introduce the correct scaling behaviour. In fact, 
if we model ^(t) as t~^+p and /i(r) as r'', then it is easy to verify that all models satisfying p -I- 2g = will give rise to 
scaling behaviour. We modify our model by allowing variations in the parameter q from the value zero. As discussed, 
the choice q = —3/2 represents a type of coherent limit. Results for such variations are shown in Fig. |2^. We see that 
as we approach the coherent hmit {q = —3/2) the value of &ioo is increased, but still doesn't come close to the data. 

We now present one final extreme variation on our standard scaling model, which does significantly improve the 
&100 problem. Fig. |2l] illustrates the effect of varying the parameter q in the limit that both the string velocity and 
the coherence length are very small. We find that in the coherent limit {q = —3/2), the matter spectrum comes close 
to the data on large scales. So in the coherent limit, the crioo normalization depends more strongly on changes in the 
parameters v and ^ than it did in the standard, incoherent case. Although this extreme model (like the mimic model 
p3| ) demonstrates that the 6ioo problem is not a necessary consequence of the assumption of scaling, we stress that our 
modifications of the degree of coherence via the parameter q as far as the coherent limit had no physical motivation. 
All known defects achieve scaling behaviour by having some kind of random decay process, such as loop production 
in the case of strings, or decay into Goldstone bosons during unwinding in the case of textures. Consequently, scaling 
seems necessarily to imply incoherence, which suggests that scaling models with constant n{t) are a special, unphysical 
limit. We see in the figure that it is only in the coherent limit that the choice of low ^ and v significantly reduces the 
6ioo problem. 



VI. MIMIC INFLATION AND VARIATIONS 

One obvious counter example to the arguments we have put forward is Turok's mimic inflation ]3^ . This is a 
scaling active model which is quite different from a standard defect model, which has been shown to produce the same 
structure of acoustic peaks in the small scale CMB as adiabatic theories. It is of interest here because it also manages 
to overcome the difficulties of the majority of models discussed in this paper by giving rise to a COBE normalized 
matter power with no deficit of power on large scales. 
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FIG. 21. Varying the scaling quality parameter q with u = 0.0 and ^ = 0.001: We plot the COBE normalized angular power 
spectrum of CMB anisotropics (left-hand graph) , and the matter power spectrum (right-hand graph) for a scaling source, with 
q = Q - long-short-dash hne, q — —1 - dotted line, q — —1.25 - short-dash line, q — —1.5 - long-dash line. Observational data 
and the prediction for standard CDM (dot-dash curve) are included for comparison. 



Mimic inflation is also of interest because the components of the source stress energy which are directly modelled 
are 9oo + © and Qd (as opposed to Qqq and 8'^ which we extract in our string model). The component Goo + 6 is 
particularly interesting because it accounts for essentially all of the contributions to the the scalar CMB fluctuations 
created at the surface of last scattering and also to the matter power spectrum. 

In this section, we briefly explore variations on the mimic model in which the ratio of Qoo + O to 0_d is allowed to 
vary freely. These variations do not respect the constraint (Ooo©"^) ^ outside the horizon. However, we believe 
that they are worth looking at because they illustrate a point about the relative dependence of ISW and matter 
contributions on Qqq + Q and Qu. 

Our mimic model is based on the version presented by Hu, Spergel and White (ref. |pG]). We take 
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We have varied the model of |Q by multiplying Od by a factor //j. The results for such variations are shown in 
Fig. ^ (for this figure we have made use of the values A = \^ Bi = 0.5, i?2 = 0.5. The long-short-dash curves show the 
standard mimic model, normalized on COBE scales. The model has an identical acoustic peak structure to standard 
CDM, though the ratio of plateau to peak height is slightly higher. The matter spectrum actually lies well above the 
observational data. If the model is normalized to CMB data on scales larger than I — 100, then matter spectrum 
lies extremely close to that of standard CDM. In other words, the mimic model matches the matter spectrum of an 
adiabatic model as well as the acoustic peak structure. 
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FIG. 22. Variations on the Mimic inflation model: We plot the COBE normalized angular power spectrum of CMB 
anisotropies (left-hand graph), and the matter power spectrum (right-hand graph) for variations on the mimic inflation model. 
The amplitude of the Qn component is multiplied by a factor fo, with fo = 1.0 (long-short-dash line), fo ~ 0.0 (dotted line) 
and fo ~ 1000 (long-dash line). Observational data and the prediction for standard CDM (dot-dash curve) are included for 
comparison. 

The short-dash curves show the results for a run where is set to be 1000 times bigger than 9oo + 0- As expected, 
this model gives rise to extremely small contributions for the matter and oscillatory components, as compared to the 
ISW. The dotted curve shows the result of a run where is set to zero. We see that the matter spectrum in this 
case is of order but slightly lower than in the basic mimic model. Therefore, the relative amplitude of matter and 
ISW contributions seems to depend in a non trivial way on the relative amplitude of and Ooo + 9, except in the 
limit that Qqq -|- 8 is small. 

In Fig. ^ we show the results of varying the Jd for different choices of the parameters A, Bi and B2, which we vary 
by multiplying each one by a factor fp from its original value. We see that with fp — 0.4, the matter power spectrum 
lacks power on large scales in the standard case of fu ~ 1.0. However, this deficit of power can be made up for by 
setting fo — 0.0. The resulting matter power spectrum is in strikingly good agreement with the observational data. 
This is in contrast to the case /p = 1.0 where the effect of reducing fu to zero is to diminish the matter spectrum. 

To summarize, the results of this section illustrate that models can be constructed which are capable of solving the 
6100 problem. However the examples of this and previous sections have also indicated that changes in the parameters 
can give contradictory results, particularly in the coherent limit. We re-emphasize that in the coherent limit there is 
no physical motivation for the concept of scaling (which comes about as a result of random processes in all known 
scaling models). Furthermore, the case which gave a good fit to the matter spectrum did not satisfy constraints on the 
superhorizon anisotropic stress that are required by isotropy. Nevertheless, these are interesting cases which should 
be investigated further to see if they could give rise to plausible models of structure formation]^ 



VII. DISCUSSION AND CONCLUSIONS 



Substantial progress has been made in understanding the predictions from cosmic defect models of structure 
formation. The modern state-of-the-art calculations |^-|ll| all represent significant advances over previous work 
1^,^-54 I^j55| , p^ , p6| , p^^ . While our 'Modeling' approach is not as explicitly tied to simulations as the work of Pen 



''As we complete this work, we learned of a preprint pll, which discusses this issue. 
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FIG. 23. Further variations on the Mimic inflation modeh We plot the COBE normalized angular power spectrum of CMB 
anisotropics (left-hand graph), and the matter power spectrum (right-hand graph) for variations on the parameters A, B\, B2- 
The amplitude of these parameters is varied from the standard case by multiplying each one by a factor /p. Also the ratio of 
6oo -I- O to Qd is varied by multiplying Od by a factor fo- We take fp = 1.0, fo = 1.0 (dot-dash line), /p = 0.4, fo = 1.0 
(dot line), /p = 1.0 fo = 0.0 (short-dash line), /p = 0.4 fo = 0.0 (short-dash line). Observational data and the prediction for 
standard CDM (dot-dash curve) are included for comparison. 



et al 1^ and Allen et al its strength is that it aUows us to probe the robustness of our key results, as a function 
of variations in the defect stress-energy histories. This lets us investigate possible variations among different defect 
models and possible systematic uncertainties in the simulations. 

We have found a serious conflict between standard scaling defect models and the current observational data. This 
conflict can be clearly expressed in terms of the '6ioo problem', where 6ioo is the bias on scales of 100/i~^Mpc. 
Current theoretical and experimental results indicate that the actual value of 6ioo is close to unity, but the standard 
defect models we considered, once COBE normalized, required &ioo ~ 5 to reconcile the predictions for the density 
field fluctuations with the observed galaxy distribution. The problem is robust over a wide variations in parameters 
describing the defect sources and over variations in fJ;, and h as well. It is tempting to exploit current uncertainties 
in the standard scaling picture to resurrect the defect picture, but we have found that truly radical deviations from 
the standard scaling assumption are required to produce a more viable model. While one variation from the standard 
model could get &ioo « 1.6 (or even of order 1 if a deviation from scaling is introduced as well), to get this close is 
done at the expense of totally failing to fit the shape of the matter spectrum and most of our ideas about how defects 
interact and evolve would also have to be wrong. 

Structure formation on 100/i~^Mpc scales is not strongly affected by changes in the matter content, and we do not 
expect that changes in the dark matter could save the scaling defect models. Probably the most interesting direction 
to take the standard defect scenario is toward an open universe or one with a non vanishing cosmological constant. 
There have been some suggestions [Q that sufficient deviations from scaling might be realized in some low VIq or 
cosmological constant models to result in a viable scenario. We also discovered a handful of very different scaling 
models (see table III) which had more reasonable values of 6ioo, but which have very little to do with the standard 
defect scenarios. They may however, represent interesting directions for future model building. 

Our results were anticipated in a number of ways by earlier work. Pen, Spergel and Turok (ref. p^) noted a serious 
bias problem for specific global defect models on scales up to 2Qh~^ Mpc and the compilation of refs. and 
in ref. looks very much like our results (and the shape of the CDM curve in ref. jl^ is almost identical to ours). 
However, in each case it was far from clear that sufficient dynamic range had been achieved to predict the COBE 
normalized density field accurately. One of the strengths of our current method is that dynamic range is not a problem 
and we feel we have covered the range of other possible uncertainties in the 'model scanning' described in sections 
III-V. In cases where comparisons are appropriate, our results are consistent with those of Pen, Seljak and Turok 
where dynamic range is also not a problem. 
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Our results are extremely negative for the standard scaling defect models and it is essential that results of this 
significance are critically scrutinized. We see four possible areas where this scrutiny might be directed. 

• Linear Einstein- Boltzmann solver : We have used standard technology to solve the linearized Einstein- 
Boltzmann equations. Initially, we modified a passive scalar code to include sources and we checked this 
with the equivalent code used in ref. The vector and tensor codes have not at this stage been checked. We 
should note, however, that an accurate Einstein-Boltzmann solver is the ultimate quantitative authority in this 
subject and that analytic fitting formulae have only limited usage if they do not agree. Calculations based on 
these fitting formulae may be wrong by up to a factor of two. 

• Modeling of two-point functions : Clearly, our results depend on our choices of defect stress-energy two-point 
functions, which were designed to fit fiat space string simulations. We have extensively probed variations about 
these two-point functions and we found that only much less motivated models are even marginally consistent with 
the data. Nonetheless, we have illustrated that some models, which at this stage are physically unmotivated, 
can fit the data on 100ft-~^Mpc. Obviously, much more work will be done on this subject. 

• Bias : It appears that 6100 ~ 5 (or even &100 ~ 3) is not consistent with the data and also that defect models 
are unable to generate such a large bias. If either of these turns out to be false the case against defects will be 
weakened. 

• Comparison with data : It is essential that we have made the comparison with the data correctly. If the current 
data turns out to be dramatically inaccurate, our conclusions would of course change. A more subtle problem 
might be the assumption of Gaussianity which goes into producing the error bars used in the comparisons 
between theory and data. We believe that the dependence of the error bars on the sample size, be it the sky 
coverage in a CMB experiment or the depth of a redshift survey, is likely to be very model dependent. Although 
the central limit theorem seems to ensure that the comparisons we make are reasonable in most defect models 
[ p8| , ^ , in more exotic cases the effect can be large. 

We have carefully scrutinized all four of these issues, and, at this stage, stand by our conclusions. However, we are 
happy to engage in critical dialogues on these topics. 

Much of the early calculations of defect models have involved various heuristic arguments, particularly concentrating 
on the effect of an individual defect on the microwave background or on matter accretion. Such approaches were able 
to give an intuitive picture of the effects of defects on structure formation, but in the modern age we feel strongly that 
arguments such as 'model X does not look so bad' must be translated into expressions for the two-point functions 
before we are willing to take them seriously — factors of two can be important and would obviously change the 
conclusions of this work. We have also demonstrated that results from coherent models can have very little to do 
with the fully active decoherent case and that it is essential that causality be strictly enforced if physically meaningful 
results are to be obtained. We are not arguing against the utility of heuristic discussions, but it is essential for the 
heuristic pictures to be brought into contact with the modern calculations if they are to have an impact. 

The complete failure of the bulk of the models we considered here is in striking contrast with the (inflation motivated) 
nearly scale invariant adiabatic models (which have no trouble getting &100 right). These two types of model are 
different in many ways and it is not necessarily easy to pin the problems on a single physical difference. Still, it is 
worth noting a number of effects which all cause problems for the defect scenarios. 

Firstly, the defect models are isocurvature models and passive isocurvature models are know to have the CDM 
perturbation at a given wavelength (relative to temperature perturbations at the same wavelength) down by a factor 
of 6 as compared with the adiabatic case. A similar effect is present in the active cases |6l), although this is not 
directly relevant to the 6100 problem, which compares temperature and CDM fiuctuations on different scales. The 
scaling (or other) assumption which links these two scales is a crucial part of the story. This difference between defect 
and adiabatic models is more directly connected with the problem of matching the galaxy data with the smaller scale 
CMB measurements, for example, within the early time window of Fig. ^, a matter which is problematic for defects, 
but which we have not emphasized here. 

Secondly, the growth of the CDM perturbations under gravitational collapse is crucial to both pictures. The fact 
that defect models lay down perturbations throughout time (as compared with the inflationary models, for which the 
perturbations are 'set up' early on) means that the average defect perturbation benefits less from the gravitational 
collapse. The temperature perturbations benefit more equally from the constant active seeding, so this is another 
physical effect which suppresses the relative perturbations in matter versus radiation. 

Finally, the previous effect is further enhanced by the generic presence of vector and tensor perturbations in the 
active models. These do not excite growing modes in the CDM, but make added contributions to the temperature 
fiuctuations, especially on COBE scales. 
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Having made these points, we must emphasize that the abihty for a causal active model to duplicate the successes of 
an adiabatic model are purely a technical matter. There is no formal 'no go theorem' against active models, it is just 
that the models which best match our current ideas about defects don't work. The best counter example is Turok's 
'inflationary mimic' model |]33|] . We have confirmed Turok's conjecture that the mimic model produces adiabatic- like 
predictions for the CDM spectrum, as well as for the CMB anisotropics (which were the original focus in [^), and we 
have added a few new scaling models to the list. We are also investigating a number of physical processes which could 
produce extremely non-scaling sources. All these models are, however, radical departures from the original standard 
scaling defect picture and they require further work to determine if they correspond to realistic physical scenarios. 

The problems with the standard defect models are likely to have a significant impact on our understanding of the 
origin of cosmic structure. The defect models were the classic examples of models of structure formation in which 
Standard Big Bang (SBB) causality holds. That is, one started with a perfectly homogeneous universe and seeded 
perturbations via causal processes in the SBB. With the demise of the standard defect models, the question arises 
whether any plausible SBB causal model exists. If the answer is negative, then this is very strong evidence for an 
inflationary origin of cosmic structure, in which SBB causality is violated to produce fluctuations of the standard 
'adiabatic' form. We are currently investigating the other possibilities of SBB causal models of structure formation 
mentioned above, but it is too early to tell if a convincing scenario can emerge. 
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